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Problem 



• Spacecraft components may be damaged due to 
airflow produced by Environmental Control 
Systems (ECS). 


• There are uncertainties and errors associated 
with using Computational Fluid Dynamics (CFD) 
to predict the flow field around a spacecraft from 
the ECS System. 

• This paper/presentation describes an approach to 
estimate the uncertainty in using CFD to predict 
the airflow speeds around an encapsulated 
spacecraft without using test data. 
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ECS System Overview 



• Environmental Control System aiaa 2002-3253 

— Prior to launch, cold air (air conditioning) flows 
downward around the spacecraft after it has been 
encapsulated in the Payload Fairing. 

— The cold air is delivered through an air-conditioning 
(AC) pipe, which intersects the fairing and flows past a 
diffuser located at the pipe/fairing interface 

- After passing over the spacecraft, it is finally 
discharged through vents 

— The Payload Fairing air conditioning is cut off at lift off. 
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ECS System Overview A 

LAUNCH SERVICES 


• Example of Environmental Control System (ECS) CFD Model 




Swirl flow 


AIAA 2002-3253 
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ECS System Overview 


• Example of an ECS system airflow test 
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Kandula, M., Hammad, K., and Schallhorn, P., "CFD Validation with 
LDVTest Data for Payload/Fairing Internal Flow," AIAA-2005-4910, 
2005. 




Overview of the Validation 

Process 
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Overview of the Validation 

Process 



• Estimate Interval within which 6 model falls with 
a given degree of confidence 

- Assume Gaussian normal distribution, 90 % confidence 


• U90%= + /-1.65*(Uval) 

Error Sources (U num, Uinput, Ud) 


u vi = ju 2 +u. f 2 -u n 2 

Val \\ num input D 


• Uncertainty Equation 


#>UCF 

Ul'ilVBRSri'V IJI CENT RAL Florida 

COLLKCF. OF KNCilMlIF-LRING, AKU COMFL TER SCIENCE 


C/™, =+/-!. 65 


90% 



JOHN F KENNEDY SPACE CENTER 


Proposed Methodology Without 

Experimental Data 


Proposed Methodology ** conservative estimate to envelop true value 


If there is no experimental data, D=0, 8 D =0, and u D =0. 
E = S — D = S 
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Ss = S-T 


E — S — D — T + fis — ( T + 8^) — S$ — Sp — Si, 


n val ^ 


num 2 "h ~h ^D 2 


- - * (f 


2 _i_ -1 1 . 2 

num ' U input 
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Report S +/-u s 

k - value (Use Student-t 
Distribution) 

Treat all input variables 
as 'random' and run 
separate CFD case 

Treat as an oscillatory 
convergence parameter 


u, 


1 


Oscillatory — o J 
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28 
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30 
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31 

30 
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41 

40 
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51 
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61 

60 
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81 

80 
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100 

1.66 

121 

120 

1.658 

infty 

infty 

1.645 



• The Rockets Behind the Missions: 

- Delta II 

- Delta IV 

- Atlas V 

- Pegasus 
— Taurus 

- Falcon 9 


• http://www.nasa.gov/centers/kennedy/launching 
rockets/ 
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EELV Public Available Information / c 

: CtWR ) 



LAUNCH SIRVICES PfiDGSAM 


Each of these vehicles have a Payload Planners Guide or Users Guide 


http://www.ulalaunch.com/site/docs/product cards/guides/DeltallPayloa 


dPlannersGuide2007.pdf 

• http://spacecraft.ssl.umd.edu/design Iib/Delta4. pi. guide. pdf 

• http://spacecraft.ssl.umd.edu/design Iib/Atlas5. pi. guide. pdf 

• http://www.orbital.com/Newslnfo/Publications/Pegasus UG.pdf 

• http://www.orbital.com/Newslnfo/Publications/taurus-user-guide.pdf 

• http://www.spacex.com/Falcon9UsersGuide 2009.pdf 
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Delta II 



• Air-conditioning is supplied to the spacecraft via an umbilical after 
the payload fairing is mated to the launch vehicle. 

• The payload air-distribution system provides air at the required 
temperature, relative humidity, and flow rate as measured 

• The air-distribution system uses a diffuser on the inlet air- 
conditioning duct at the fairing interface. 

• If required, a deflector can be installed on the inlet to direct the 
airflow away from sensitive spacecraft components 

• The air can be supplied to the payload between a rate of 

• 1300 to 1700 scfm. 

• Diameter of Fairing is 3meters 
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Delta II - Continued 


Fairing Wall 
(Inside} 


Air Flow 


Ai ^Conditioning , 
Inlet Diffuser 


Air Flow 
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HHMWIIREUD 

Lanya'd 

Disconnects 


Air- Conditioning 
Duct 



Air-conditioning duct and diffuser 
system is ejected at liftoff 


Figure 4-1, Payload Air Distribution System 
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AirOAuftOf 


Acous:*: Blankets 



Fsirlng Composite Shefl 


Air-Conditioning Ouct 


Figure 4-1. Standard 4-m Composite Fairing and S-m Composite Fairing Air-Conditioning 

Duct Inlet Configuration 

Air flows around the payload 
and is discharged through 
vents in the aft end of the 
fairing. 

Fairing sizes 4meter and 5 
meters in diameter 
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The air is supplied to the payload at a 
maximum flow rate of 36.3 kg/min to 72.6 
kg/min (80 to 160 Ib/min) for 4-m fairing 
launch vehicles and 90.7 kg/min to 136.0 
kg/min (200 to 300 Ib/min) for 5-m fairing 
launch vehicles. 



LV273* 

<tSop 


View Q 
Not to Scale 

Figure 4-2. 5-m Metallic Fairing Payload Air-Distribution System 
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• Internal ducting defectors in the PLF direct the gas upward to prevent direct 
impingement on the spacecraft. 

• The conditioning gas is vented to the atmosphere through one-way flapper doors 
below the spacecraft. 

• The PLF air distribution system will provide a maximum air flow velocity in all 
directions of no more than 9.75 mps (32 fps) for the Atlas V 400 and 10.67 mps (35 
fps) for the Atlas V 500. 

• There will be localized areas of higher flow velocity at, near, or associated with the 
air conditioning outlet. 

• Maximum air flow velocities correspond to maximum inlet mass flow rates. 

• Reduced flow velocities are achievable using lower inlet mass flow rates. 

• Flow Rates 

A) Atlas V 400: 0.38-1.21 kg/s ±0.038 kg/s (50-160 Ib/min ±5 Ib/min), 

B) Atlas V 500: 0.38-2.27 kg/s ±0.095 kg/s (50-300 Ib/min ±12.5 Ib/min) 

• Fairing sizes are 4meters and 5 meters in diameter 
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Pegasus 



• The fairing is continuously purged with filtered 
air. 


• The flowrate of air through the fairing is 
maintained between 50 and 200 cfm. 


• The air flow enters the fairing forward of the 
payload and exits aft of the payload. There are 
baffles on the inlet that minimize the 
impingement velocity of the air on the payload. 

• Fairing diameter is 0.97 meters 
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Taurus 



• Upon encapsulation within the fairing and for 
the remainder of ground operations, the 
payload environment will be maintained by 
the Taurus Environmental Control System 
(ECS). 

• Fairing inlet conditions are selected by the 
Customer 


• Fairing diameters are 63 inches and 92 inches 
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Falcon 9 



• Once fully encapsulated and horizontal, the 
Environmental Control System (ECS) is 
connected 


• Payload environments during various 
processing phases are: 

- In hanger, encapsulated - Flow Rate: 1,000 cfm 
— During rollout: 1,000 cfm 

- On pad: Variable from 1000 to 4500 cfm 


• Fairing diameter is 5.2 meters 
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(3) Configurations 



• Fairing Sizes are approximately lm, 1.6m, 2.3m, 3m, 
4m, 5m in diameter. 

• (3) generic fairing diameters are selected to envelop 
the EELV fairing configurations 

- 0.75m 

- 3.5 m 

- 5.5 m 

• Inlet Conditions range from 1000 cfm to 4500 cfm 

• Spacecraft diameters range with fairing sizes, a generic 
spacecraft was drawn and scaled accordingly 
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(3) Configurations 



• CAD model of the spacecraft was created in 
Pro/ENGINEER, 0.75m 
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(3) Configurations 
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• Snappy Hex- Mesher 


0.75m 
Configurati 
(6762865 nurr 
of cells) 
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3.5m Configuration 
(8594480 number of 
cells) 


5.5m Configuration 
(6980673number of 
cells) 


COLLKCF. m : ICNCilNREUtl^i AN'U CDSEFL TEH SCIENCE 



JOHN F KfNNFDY SPACE CENTER 


• OpenFoam 



0.75m 

Configurati )n 
(6762865 number 
of cells) 
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CFD Modeling 



SimpleFoam 



on©_coarsegrid (m/s) 



3.5m Configuration 
(8594480 number of 
cells) 


5.5m Configuration 
(6980673number of 
cells) 
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Uncertainty Calculation 
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Proposed Methodology 


S u val 
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num ■ u 


input 


InputH/ariableE 

Descriptions 

BiasE [ 

Grids 

BISridsIlonsideredllordachllonfigurationS 

sh 

InletS/elocityS 

BoundarytEonditionSowdndUiighS 

10%S j 

OutletGIPressureS 

BoundarylEonditionlowllndliighS 

2%s ; 

Turbulenced/lodelS 

SA,de-realizable,[lcwSSTS 


WallGFunctionsS 

withdndSvithoutS 

sh 

RoughSA/allS 

Functions 

smoothes. SoughS 

m 

CompressibilityS 

incompressiblels.lompressibleS 

50 

SolverS 

OpenFoam,[¥luent,[MARCCM+S 

50 

FluidGiPropertiesS 

kinematicl^iscosityliulepresents[lir[|0-50-100][SlegS 

CS 

1.36,1.5,2.306e- 
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Uncertainty Calculation / 


LAUNCH StRVlCES 


Case# 



Parameter 

Grid 


u vai = 1.746 * 

UCF 
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0.75 3.5 5.5 


1 
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med 
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3 

Boundary Conditions 
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inlet velocity low 

4 

4 
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inlet velocity high 
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5 

5 

6 
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6 

6 

7 


pressure outlet high 

7 

7 

7 

Turbulence Models 
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8 

8 

a 
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9 

9 
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10 
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Wall Functions 

without wall functions 
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12 

12 
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13 
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15 

15 

Fluid Properties 
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16 
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Results 0.75m Configuration J c 

MTER ) 
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Solution and Uncertainty Contour Plots 
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Results 0.75m Configuration J c 
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Solution Line Plot 
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• Uncertainty Line Plot 
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Results 0.75m Configuration J c 

MTER ) 
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Uncertainty Ranking 


L 



• The uncertainty for each of the input variables were ranked by the non- 
dimensionalizing the difference in the results by the freestream value and 
ranking from greatest uncertainty to least uncertainty. 
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Input0/ariable0 

Description]?] 

Bias0 

Mean!?] 

Velocity]?] 

Uncertainty!? 

(m/s)0 

MeanSslon- 

DimensionalizedC? 

Uncertainty0 

Normalized5 

Ranking5K>0 

Numbered! 

Ranking0 

Grid0 

3Cgrids0 

considered!?] 

50 

1.62870 

0.05430 

13.400 

20 

Inlet0/elocity0 

Boundary!?] 

Condition!?] 

io%0 

1.31150 

0.047370 

11.690 

50 

Outlets 

Pressure!?] 

Boundary!?] 

Condition!?] 

2%m 

1.14780 

0.03830 

9.450 

80 

Turbulence!?] 

Model!?] 

SA,Qke- 

realizable,0 

kwSST0 

50 

1.46280 

0.04880 

12.040 

40 ; 

\A/all[3Functions0 

with!]&nd0 

without!?] 

m 

0.82860 

0.02760 

6.810 

90 

Rough0/Vall0 

Function!?] 

smoothlsJ 

rough!?] 

50 

1.52370 

0.05080 

12.530 

30 

Compressibility!? 

incompressible!? 

vs.0 

compressible!?] 

50 

1.31280 

0.04380 

10.810 

60 

Solver!?] 

OpenFoam ; 0 

Fluent,!?] 

STARCCM+0 

50 

1.6730 

0.05580 

13.770 

10 

Fluid!?] 

Properties!?] 

kinematic!?] 
viscosity!lhu0 
represents!3&ir0 
[0-50- 100] SI eg!?] 

cm 

1.36,1. 5,2. 306e- 
050 

1.15360 

0.03850 

9.500 

70 ! 
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Results 3.5m Configuration 
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Solution and Uncertainty Contour Plots 
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one_coarsegrid (m/s) 
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Results 3.5m Configuration 


Solution Line Plot 
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Uncertainty Line Plot 
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Uncertainty Ranking 


Configuration 
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Input Variable 

Description 

Bias 

Mean 

Velocity 

Uncertainty 

(m/s) 

Mean Non- 
Dimensionalized 
Uncertainty 

Normalized 
Ranking % 

Numbered 

Ranking 

Grid 

3 grids 
considered 


0.6829 

0.0228 

8.28 

7 

Inlet Velocity 

Boundary 

Condition 

10% 

0.7919 

0.0264 

9.59 

6 

Outlet 

Pressure 

Boundary 

Condition 

2% 

1.4606 

0.0487 

17.70 

1 

Turbulence 

Model 

SA, ke- 
realizable, 
kwSST 


1.3487 

0.045 

16.35 

2 

Wall Functions 

with and 
without 


0.6139 

0.0205 

7.45 

9 

Rough Wall 
Function 

smooth vs. 
rough 


1.0531 

0.0351 

12.75 

3 

Compressibility 

incompressible 

vs. 

compressible 


0.8252 

0.0275 

9.99 

5 

Solver 

OpenFoam, 

Fluent, 

STARCCM+ 


0.841 

0.028 

10.17 

4 

Fluid 

Properties 

kinematic 
viscosity nu 
represents air 
[0-50-100] deg 
C 

1.36,1. 5, 2.306e- 
05 

0.6345 

0.0212 

7.70 
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Results 5.5m Configuration 
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Solution and Uncertainty Contour Plots 
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Solution Line Plot 
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• Uncertainty Line Plot 
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Uncertainty Ranking 


Configuration 
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Input Variable 

Description 

Bias 

Mean 

Velocity 

Uncertainty 

(m/s) 

Mean Non- 
Dimensionalized 
Uncertainty 

Normalized 
Ranking % 

Numbered 

Ranking 

Grid 

3 grids 
considered 


2.0203 

0.0673 

12.44 

3 

Inlet Velocity 

Boundary 

Condition 

10% 

1.6198 

0.054 

9.98 

6 

Outlet 

Pressure 

Boundary 

Condition 

2% 

2.0173 

0.0672 

12.42 

4 

Turbulence 

Model 

SA, ke- 
realizable, 
kwSST 


2.3049 

0.0768 

14.19 

i 

Wall Functions 

with and 
without 


1.4902 

0.0497 

9.18 

7 

Rough Wall 
Function 

smooth vs. 
rough 


1.4901 

0.0497 

9.18 

8 

Compressibility 

incompressible 

vs. 

compressible 


1.4256 

0.0475 

8.78 

9 

Solver 

OpenFoam, 

Fluent, 

STARCCM+ 


1.8172 

0.0606 

11.20 

5 

Fluid 

Properties 

kinematic 
viscosity nu 
represents air 
[0-50-100] deg 
C 

1.36,1. 5, 2.306e- 
05 

2.05 

0.0683 

12.62 

2 
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Comparison to Previous LDV Test j c 

R * ^ 
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Kandula, M., Hammad, K v 
and Schallhorn, P v "CFD 
Validation with LDV Test Data 
for Payload/Fairing Internal 
Flow/' AIAA-2005-4910, 2005. 


Variable 

Bias 

Velocity Inlet 

3% V 

Kinematic Viscosity 
[0-100] Deg C 

[1.36, 1.50, 2.306] e-5 m2/s 

Pressure Outlet 

3% 

Turbulence 

ke-realiazable, kwsst, SA 
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Velocity (mA) 6 Velocity (m/») 
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Comparison to Previous LDV Test 
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Ut (x=9in) 

*3 
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— Uncert_low 
— Uocert_Hifh 
-♦-CxpenmenUl 
— CfO 


OS 



IS 


i (meters) 

Uy(t*-Sin) 



Assumes 99% Confidence 
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Comparison to Previous LDV Test J c /- 

? Vy 
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Uz (x=- 3in) 


x (meters) 
Uz (z=*2in) 


2.5 



— Oncert_low ^ 

— Unceft_Hifh > 0 5 
“•“CxpenmcntJ | 

06 -C« * fl 

-05 


*15 


Uy (z=-3in) 



a (meters) 

Uy(z=-2in) 
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— Oncerr_H*h 
“•"£*penmef»Ul 
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Comparison to Previous LDV Test j c 
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Assumes 99% Confidence 
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Conclusion / Recommendation A 
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Proper validation with experimental data should oe 
used to verify ECS impingement requirements 



• This research proposes a CFD uncertainty 

methodology when experimental data is unavailable 
and unobtainable 


• Couples Student-T Distribution to the number of CFD 
models and input parameters 

• All input parameters considered had the same order of 
magnitude uncertainty 
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